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Abstract 

We extend the matrix decomposition method(MDM) in classifying the 2xN xN 
truly entangled states to 2 x M x system under the condition of stochastic local 
operations and classical communication(SLOCC). It is found that the MDM is quite 
practical and convenient in operation for the asymmetrical tripartite states, and an 
explicit example of the classification of 2 x 6 x 7 quantum system is presented. 



1 Introduction 

Entanglement is an essential feature of quantum theory, describing a quantum cor- 
relation that exhibits nonlocal properties. In the seminal work [1], Einstein, Podolsky, 
and Rosen (EPR) demonstrated through a Gedanken experiment that the quantum me- 
chanics (QM) can not provide a complete description of the "physical reality" for two 
spatially separated but quantum mechanically correlated particles state which is now 
known as entangled state. The subsequent Bell theorem manifest the nonlocal character 
of the quantum correlation in the violation of Bell's inequalities [2]. As the quantum 
information science develops, the impact of entanglement goes far beyond the testing of 
the conceptual foundations of QM. Entanglement is now of central importance in the 
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quantum information theory (QIT) and is thought as the key physical resource to real- 
ize quantum information tasks, such as quantum cryptography [31 H] , superdense coding 
[5l [6], and quantum computation etc. This necessitates the qualitative and quanti- 
tative description of the entanglement |8]. However due to the lack of suitable tools for 
characterizing the entanglement, very limited quantum state space was explored in the 
quantum information theory. 

In quantum information processing (QIP), two states are suited to implement the 
same task if they can be mutually converted by Stochastic Local Operations and Classi- 
cal Communication (SLOCC) [9], and therefore they are said to be in the same equivalent 
class. For three qubits, known result is that there are two kinds of true tripartite en- 
tanglement classes for pure state, namely, GHZ and W states [9j. As the dimensions of 
each party increases nontrivial aspect shows up, i.e., non-local parameters may resides 
in the entangled states of 2 x x system when A^ > 4 [TOl [11] . Many investigations 
concerned the classifications of 2 x M x A^ states has been done in [TUl [121 [IS]- In the 
Refs. [ini [12] , an iterated method was introduced to determine all the inequivalent classes 
of the entangled states of 2 x M x A^ system based on the "range criterion", where the 
entanglement classification of the low dimension system is a prerequisite for the high di- 
mensions ones. Practical classifications of dimensions up to 2 x 4 x 4 and the related 
systems of 2 x (M + 4) x (2M + 4) were given in [lOj . With the increasing of dimensions, 
the complexity of the method grows dramatically because of the iterated nature of their 
inequivalent proof of the entanglement classes. In a recent work [TT] a novel method of 
classifying the pure state of 2 x N x N systems was introduced in which all the inequiv- 
alent true tripartite entanglement classes can be determined directly by using merely the 
elementary operations on the cubic grid form of the state. 

The present work deals with the more general case: quantum state of 2 x M x A^ 
systems (pure state if not specified). We show that the method we introduced in [TT] 
can be generalized to the classification of true entangled states of 2 x M x A^ systems. 
And therefore all the inequivalent classes can be generated directly and no followed-up 
inequivalence proof of these classes is needed. The content goes as follows, in section 2, 
by representing the 2 x M x N state in the form of matrix pairs, the 2 x M x N states are 
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Figure 1: The cubic form for 2 x 3 x 4 state. 



divide into inequivalent sets under SLOCC. The detailed classification procedures with 
these inequivalent sets are presented in section 3 and a concrete example of classification 
of 2 X 6 X 7 system is given. Finally, in section 4 we give some concluding remarks. 

2 Matrix pair representation of 2 x M x state 

Adopt the conventions of [iTj, an arbitrary state of 2 x M x can be written as 

|^2xMxAr) = Tr[r{ij- fc} ll)2®'^\® V'J] , (1) 

where, represents the first qubit, ipi and ?/'2 has the dimension of M and N separately; 
r{ij,A;} and r{2j',fc} are M x complex matrices (we assume M < N without loss of 
generalities). Then the state can be written in the following compact form 



|^2xMxAf) 



r Ti 1 

T2 



(2) 



Clearly, to every state of 2 x M x A^, there is a form of Eq.([2]) that corresponds to it, and 
a pictorial description of the state is straightforward, see Fig.([T]). 

The reduced density matrix of state ^2xMxn is defined as pa = Tr^^|\E')(\E'|, where 
a, /3, 7 can be '?/'o, V^i, '^2- For three-partite systems, ture (genuine [9j) entanglement means 
that the determinant of the reduced density matrix of each partite is nonzero. For the 
2 X M X N systems, this is equivalent to r^p^^) = 2, r{p^^) = M, r{p^^) = N. The density 
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matrix in the form of the matrix pairs can be expressed as 

where i, i' = 1, 2; j, j' = 1,2, ■ ■ ■ , M; k,k' = 1, 2, • • • ,N. Then the reduced density matrix 
(take ip2 as an example) is 

= ^i^i)jk'{^i)jk 
ij 

= Erlr.- (4) 

i 

We can infer that if r{p^^) < N (dct(p,02) = 0), there will be ILOs that transform p^^ to 
who has at least one column or one row of zeros. Without loss of generalities suppose 
the kth column of p'^^ are zeros, for the element {k,k) of p'^^ we would get 

ij 

= 0, (5) 

which indicates that (r^)^^ = for all i and j. In the pictorial description of the 2xM xN 
state, this corresponds to the case that the cubic grid has a whole plane of zero coefficients. 
Clearly in this condition the entanglement oi 2 x M x N system reduces to the case of 

2 X M' X N' with M' < M or/and N' < N which should in principle be considered as an 
entanglement system of 2 x M' x N'. 

3 Classification of 2 x M x TV State 

Two 2xM xN states 4f and * are said to be SLOCC equivalent if they are connected 
via invertible local operator (ILO) . That is ^ is SLOCC equivalent to ^ if 

|*2xMxAr) = P (g) (5 |^2xMxAr) , (6) 

where T, P, Q are invertible complex matrices of dimension 2x2, MxM, and NxN which 

act on ipQ, ipi, "4^2, respectively. Neglecting the extra factor of the determinant of matrices. 
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T, P, and Q correspond to the special linear groups of SL{2, C), SL{M, C), SL{N, C) [9]. 
Take |\E'2xMxAf) in the form of Eq.(l2]), the ILO operators T, P, Q in Eq.(l6]) can be written 



as 



1^ 



2xMxN 



tu ti2 




' PTiQ ' 


^21 ^22 




[ PT^Q J 



(7) 



From Eq.(j2]) and Eq.(I7j) we can see that the SLOCC equivalence of the quantum state turns 
to the connectivity of the matrix pairs (Fi,F2) under the special linear transformations 
T,P,Q. Define the set that contains all the matrices pair (Fi,F2) as C. The whole space 
of C can be partitioned into numbers of inequivalent sets with different n, I 



Cn, 1 = {{Ti, F2)| rmax(aiFi + /31F2) = n, rmin(a2Fi + /32F2) = /} 



(8) 



where r^ax and rmin represent the the maximum and minimum rank of the matrices 
respectively, and we let r denote the rank of matrix hereafter; aj, (3i E C and \ai\ + \f3i\ 7^ 0; 
I e [0,n], n e [0,M]. 

3.1 Classification on sets C„ ^ with n = M 

We start our classification of Cnj in2x M x N system from the case n = M. Our aim 
is to construct the subsets cm,i which: (i), it includes representative elements of all the 
inequivalent entanglement classes; (ii), each inequivalent class has only one representative 
element in cm,i- 

V (Fi, F2) G Cm, I there always exists an ILO operator T that 

f Fi 1 _ C til ^12 1 f ^1 1 /„x 
I ^2 j t21 t22 ) I ) 

where r(tiiFi + ti2F2) = M, r(t2iFi + t22F2) = /. So we just assume r(Fi) = M and 
'"(F2) = / without loss of generality. Two specific ILOs P and Q can transform (Fi,F2) 
into the following form 

Fi \ P,Q, 



( Emxm 'Omx{N~M) ) 
( AmxM Bmx{N-M) ) 



(10) 



where E is an unitary submatrix, hereafter represents zero submatrix, A and B have 
the same dimensions as E, separately, and all of them have the subscripts as their 
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dimensions. If (A^ — M) > M, even if its submatrix B has the maximum rank of M, the 
right hand of Eq. flTO|) can still be further transformed by ILOs into 



r2 



EmxM 



JMx{N-2M) ^MxM 
EmxM 



(11) 



JMx{N-2M) 

In the form of the cubic grid (Fig.([T])), this corresponds to that {N — 2M) vertical planes in 
the middle of the cube are zero planes, which is actually an entangled states of 2 x M x 2M 
according to the statement below Eq.Q. Thus here we only consider the case M > N/2. 

For arbitrary matrix pair with the form of the right hand of Eq. ffTOj) . we can implement 
the following transformation by ILOs 
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where r(£") = r{B) and -E" is submatrix with typical form of {Sir{E')x{{n-m)-r{E')),Er[E')xr(E'))i 
A', EiA' are also square submatrices with the dimensions (m — r{E')) x [m — r{E'))] 
the rest of the matrices are partitioned accordingly, i.e., Oi^', B' have the dimension 
(m — r{E')) X r{E'), Oia, 02a have the dimension of (m — r{E')) x [n — m), Oih, O26 have 
the dimension r{E') x (m — r{E')), Ei', O2C have the dimension r{E') x r{E'). After this 
transformation Fi = {Emxn, 'Omx(n~m)) is unchanged, F2 becomes a quasidiagonal matrix 
and we named this procedure step i. 



Next we repartitioned the matrices on the left hand side of Eq.f lT2|) as follows 
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This is named as step ii. Consider the submatrix B', if it is not identically zero we can 
perform the transformation of step i on the left top submatrices of Eq. f|T3|) 

/ 



( EiA' Ola ) 

{A' B' ) 



step i 
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This procedure can be done repeatedly (suppose repeat n times), until the B^"'^ is iden- 
tically zero or has zero dimension. We can get that the matrix pair (ri,r2) can be 
transformed into the following form 



( 
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where the transformed Fi is just (-Emxa/, OA/x(7v-Af))5 and Ei^2 are defined according to 
the partition lines. 

As a concrete example here we show how this whole procedure is proceeded on the 
sets of 6*4,; of 2 X 4 X 6 state. The transformation of Eq. (fT2l) is start with 
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where 
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Here the rank of -84x2 must be 2 due to the argument below Eq. ffTT]) . The step ii goes as 
follows 
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Next we repeat the step i to the up-left submatrices of the right hand side of Eq. (lT9|) . 
This iteration of step i depends on the rank of B' . 



(1), r{B') = 0. In this case the matrix pair (ri,r2) become 
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r2 
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And there are three different forms, i.e., 
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(1.1) 
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,(1-2) 
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correspond to two Jordan canonical forms of A', J 



,(1.3) 
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(2), r{B') = 1. In this case 
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where A" , B" are matrices of 1 x 1 and E" = (0, 1). Again apply step i on {A" B") we 
have 



(2.1), r{B") = 
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(2.2) r{B") 
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For Eg. (12^ . A" is equivalent to the case of A" = according to theorem 1 of [TT]. For 
Eq. fP^l) . in the next step of step ii, B^^'' will be a matrix of dimension zero, and satisfies 
r(B^^^) = 0, thus the procedure is stopped. We get two inequivalent forms of r2 
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(3). r{B') = 2. In this case 
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Thus here is only one class, where r2 has just the form of Eq.( l27I) . In the following, we 
shall see that these six cases correspond to the six inequivalent entanglement classes in 
2x4x6 systems, which agrees with the result of Ref . [I2] . 



In all, for every (Fi, F2) G Cm, i, there exists an ILO transformation that make 
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maps Cm, I to cm,i, where cm, i ^ Cm, i and 
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has the form of Eq.(^. Eq.([28 
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(29) 



) El ) ' ^ \^ E2 
Thus we have separated the classification of Cm,i into two procedures: 1, the construction 
of E2 matrix; 2, classification of J. And for the second procedure, we have already 
complete the classification in pjj. We have the following theorem 

Theorem 1 V (Fi,F2) G cm,i, the set cm,i is of the classification of Cm,i- (i) if two 
states are SLOCC equivalent then they can be transformed into the same matrix vector 
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(ri,r2); (ii) this matrix vector is unique in Cm,i, that is z/(ri,r2) is SLOCC equivalent 
with (Fi, r2), then (Fi, Fg) = (Fi, F2), (F2 = F2 means that E2 = E'2 and their Jordan forms 
of J are equivalent under the condition of theorem 1 Ref.fTl^ ) 

Proof: 



(i) The proof of this statement is straightforward, since in every step of transformation 
only invertible operators take part in. 



(ii) Suppose 
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First we show that the T' transformations can always be replaced by ILO operators 
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Because F2 in (Fi,F2) G cm,i has a form of direct sum of J and E2 as shown in the 
definition (l29i) . Thus when the dimension of J does not equal zero, there are no zeroes in 
pivot of T' and the left hand side of Eq. fl3T]) can be separated into two parts 
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is the LU decomposition of T' [\A\ ; Ei has the same dimen- 



For the J sub-matrix we have proved fTT] there exists Pj, Qj which make 
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where A' = ^. It is simple to verify that such kind of Px,y, Qx,y satisfying the equations does 
exist (see Appendixes of [TT] for detailed derivations). Thus Pc = PxPy and Qc = QyQx 
will make 



1 
A 1 



a (3 
7 



PcEiQc 
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Combine Eq.f lMI) and Eq.( l36i) we can get such Pq = Pj Q) Pq, Qq = Qj © Qc that satisfy 
the following equation 
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which is just Eq. (!3Tl) . 

In one word, Eq. fj30|) can always be written as 
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The effective transformation P",Q" that keep Fi invariant must be of the form 
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where Det(y) 7^ 0. The transformation of F2 reads 

P'TaQ" = P"{A, B)Q" = {P"AP"-^ + P"BX, P" BY) 
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(40) 



where A is the M x M submatrix, and B is the M x (A^ — M) submatrix. Since P" and 
Y both are ILO operators, the rank of submatrix B, is unchanged and it can be further 
transformed to form of Eq. (fT2l) 

A' B' 



E' 



(41) 



We get that if two states are SLOCC equivalent then E' block of Fg and F2 must be 
identical. In Eq. (HT|) we see that E' block is diagonalized in F2 which guarantee that we 
can proceed our proof to E" ( E^^\ E^^^ and so on) by leaving the E' block unchanged. 
The rest steps of this proof are just like argument after Eq.l HOj) . and finally we can get 
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that if (Fi, Fg) is SLOCC equivalent with (Fi, F2) then Fg and F2 have the same canonical 
form in the set of Eq. fl29l) . (However there exists the special case that the dimension of J 
equals zero, in this case there can be zero elements in the pivot of the nonsingular square 
matrix T' . T' can then be decomposed as decomposed as 
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where a, /3, 7, A G C, Pt' = 
equation are nonsingular. It can be s 
Pz, Qz which act on Fi and F2, i.e. 



El 
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and both matrices on the righthand side of above 
low that Pt' can be compensated by some operators 
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see Appendixes of ^Ij) Q.E.D. 

3.2 Classification on sets C„ ^ with n = M — 1 

For the set of cm-i,« G Cm-i,i, the construction of cm-i,i is essentially the same as 
that cn-1,1 of 2 X N X N systems in [TT]. Here for the sake of convenience we only take 
a 2 X 7 X 8 state as a demonstration, i.e., Cm-i,i = Ce,/- The matrix pair (Fi,F2) can be 
transformed into the following form 
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where r2 can then be expressed as 
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Further simphfication can be proceeded according to the vector(or submatrices) c, r. 
There are four cases in general, i.e., (1), (c = 0,r = 0); (2), (c 7^ 0, r = 0); (3), 
(c = 0,r 7^ 0); (4), (c 7^ 0,r 7^ 0). Here c 7^ means that r(c) > 1 and different 
ranks will result in different classes, i.e.. 
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Clearly, analogous with the set Cn-i,i in Ref. [H], we can finally get the following set 

CM^i,, = {(A,r)|r(r) = /;r= Q iA,T)eCM-i,i}. (49) 

J represents the Jordan canonical form. 
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Theorem 2 V (A, F) G C]\i-ij, the set cm-i,i is of the classification of Cm-i^i. (i) 
suppose two states are SLOCC equivalent, they can be transformed into the same matrix 
vector (A,r); (ii) this matrix vector is unique in cm-i,i, that is suppose (A,r') is SLOCC 
equivalent with (A,r), (A,r') = (A, F) (F' = F means Js are equivalent in the condition 
of theorem 1 and B' = B). 

The proof of this theorem is in the same manner as that of theorem 2 in [IT] . 

Following the procedure of n = (M — 1) and the method introduced in [11], we can 
further classify the more generally cases, i.e., n = M — i with i > 1. Here we neglect 
the proof and only give a concrete example of the classification of this kind. We give a 
complete classification of 2 x (M + 5) x (2M + 5) for M = 1, i.e., 2x6x7 state whose 
classification has not been presented in literature so far. 

Classes of sets Cg,; : for all inequivalent classes in cg,/, they have the same form of Fi in 
the definition ^ 
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1 

1 





1 

1 

1 

1 

1 






1 
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1 























1 
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Here the square matrices of {x}„xn in Eq.(l5 H [52]) consists of all the inequivalent classes 
of sets c„,z in 2 X n X n states. For example the first matrix of Eq.f lSTjl is made up by all 
the genuine entanglement classes of the sets c^^i in 2 x 5 x 5 state and plus the one with 
{xjsxs = 0, thus there are (26 + 1) [15] inequivalent forms of this matrix. 

Classes of set c^^: for all inequivalent classes in c^j, they has the same form of A 

" 1 " 

1 

1 

1 

1 





A 



The different r2S are 

X X 

X X 



1 

1 

1 







1 

1 

1 

1 





1 



1 

1 

1 

0' 

1 

1 
1 
1 

1 _ 

" 

1 
1 
1 
1 
1 



1 

1 



1 

1 

1 

1 0' 

1 



1 

1 

1 



(53) 



(54) 



(55) 



(56) 



Same as that of cqj, {x}2x2 here has three different forms. 

Thus we get (26 + l) + (13 + l) + (5 + l) + (2 + l) + l + l + (2 + l) + l + l + l + l + l + l = 61 
inequivalent entanglement classes in 2 x 6 x 7. It is clearly to see that this method is 
simple and effective, meanwhile each entangled state can be read out directly from the 
matrix pairs. 
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4 Conclusions 



In summary, we have generalized our method of entanglement classification under 
stochastic local operations and classical communications to the more general case of 2 x 
M X N systems. Two examples of 2 x 4 x 6 and 2x6x7 are given where all their 
inequivalent entanglement classes are determined. Because the classification procedure is 
essentially a constructive algorithm, the method can serve as a powerful tool in practical 
entanglement classifications with the aid of computers. Most importantly a wide range of 
state space is explored which provide a rich resource for possible new applications in the 
quantum information theory. 
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